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Q ! ABSTRACT 

o 

r^i ', I consider analytically particle acceleration at relativistic shocks, in the limit of pitch 

pP angle diffusion and large shock Lorentz factors. The derived energy spectral index (k = 

(1 + \/l3)/2 ~ 2.30), and particle pitch angle distribution at the shock are successfully 
compared with the results of numerical solutions. This totally analytic derivation is 
completely independent of the detailed dependence of the diffusion coefficient D(fj,,p) 
J> ■ on either parameter, and it is argued on a physical basis that the orientation of the 

magnetic field is also irrelevant, making k a universal index. 
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■ 1. Introduction 

Q-f 

There is currently a growing interest in the acceleration of non-thermal particles at highly 
relativistic shocks. There are three classes of bona fide relativistic sources: beyond the well- 
established extra-Galactic (Blazars) and Galactic (superluminal) sources, both of which exhibit 
superluminal motions, it is now also well-established that Gamma Ray Bursts (GRBs) display 
highly relativistic expansions, with Lorentz factors well in excess of 100. Other classes of relativistic 
sources may include Soft Gamma Ray Repeaters (SGRs), whose recurrent explosions are largely 
super-Eddington, and special SNe similar to SN 1998bw, which displayed marginally Newtonian 
expansion (~ 6 x 10 4 km s _1 ) when optical emission lines became detectable, about a month after 
the explosion. 

With the discovery of GRBs' afterglows, it has now become feasible to derive the energy spectral 
index k of electrons accelerated at the forward shock, as a function of the varying (decreasing) 
shock Lorentz factor 7, provided simultaneous wide-band spectral coverage is available. With the 
launch of the USA/Italy/UK mission SWIFT, these data will become available for a statistically 
significant number of bursts, testing directly models for particle acceleration at relativistic shocks. 
Furthermore, since GRBs must also clearly accelerate protons, the same index k may determine 
the spectrum of ultra high energy cosmic rays observed at Earth. 



However, until recently, both the lack of astrophysical motivation and the difficulty inherent 
in treating highly anisotropic distribution functions have stiffened research on this topic. Early 
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work, both semi-analytic and outright numerical, has concentrated on barely relativistic flows with 
Lorentz factors of a few, well suited to Blazars and Galactic superluminals, but clearly insufficent 
for GRBs, the only exception being the numerical simulations of Bednarz and Ostrowski (1998). It 
is the purpose of this Letter to perform an analytic investigation of the large 7 limit, to establish 
which (if any) of the properties of the particles' distribution function depend upon the physical 
conditions of the fluid. 



I deal first with pure pitch angle scattering, and then (Subsection 2.3) I discuss oblique shocks. 
In the well-known equation for the particles' distribution function, under the assumption of pure 
pitch angle scattering, 



/ is computed in the shock frame, in which are also defined the distance from the shock z, the fluid 
speed in units of c, u, and fluid Lorentz factor 7. Instead, the scattering coefficient D, particle 
momentum p and particle pitch angle cosine, fi, are all defined in the local fluid frame. I make no 
hypothesis whatsoever about D, except that it is positive definite and smooth. We place ourselves 
in the shock frame, and call z = the shock position; the upstream section is for z < 0, so that the 
fluid speeds are both > 0. The above equation admits of an integral: by integrating over fi and z 
we see that 



independent of z. The required boundary condition for /, i.e., that / — ► as z — ► —00, implies 
that the constant, upstream, is 0. Downstream, Eq. 2 is also a constant, but, because of Taub's 
jump conditions, it is not the same constant as upstream. Since the required boundary condition 
for / far downstream (foo) is that it becomes isotropic and / is a relativistic invariant, we see that, 
far downstream, nfocd^ = 0. We thus have 



2. The analysis 




(1) 




(2) 




z < 
z > 



(3) 



where the inequality in Eq. 3b (which will become necessary later on) derives from the obvious 
constraint / > 0. 
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2.1. Upstream 



/: 



I begin the analysis by considering Eq. 3a in the limit of very large shock Lorentz factors, in 
which case, upstream, u — > 1. For u = 1, this reduces to 

\i + n)f4i = o. (4) 
1-1 

Since 1 + fi > everywhere in the integration interval except of course at fi = —1, and since 
/ > we see that, for u = 1 we must have / oc <5(/z + 1) where 5{x) is Dirac's delta. Thus, in 
this limit, the angular dependence factors out. For reasons to be explained in the next Subsection, 
we shall also need / for 1 — u <C 1, but still ^ 0. To search for such a solution, we let ourselves 
be guided by the solution at u = 1: thus, we let the angular dependence factor out, and use the 
Ansatz f = g(z)w((fj, + l)/h(u)). Here h(u) is an as yet undetermined function of the pre-shock 
fluid speed such that h(u) — ► as u — > 1. In this way, as the speed grows, the angular dependence 
becomes more and more concentrated toward fj, = — 1, as required by the previously found solution 
for u = 1. Introducing our Ansatz into Eq. 1 I find 

7 dg _ 1 d / ^ ^i 2 )^^ — ^-D-i^ 2 /qn 

where I defined D_i = D(fi = — l,p), and I factored the eigenvalue for future convenience. Con- 
centrating on the angular part, and defining (fj, + l)/h(u) = y, and w = dw/dy, w = d 2 w/dy 2 , I 
find 

— ^ — " + W)^ ( (Ai,p)( " H) — w — = (6) 

We are interested in a solution of the above equation only in the limit h(u) — > 0, the only one in 
which our factored Ansatz is a good approximation to the true /. In this case, the term oc w is clearly 
subdominant, and can be neglected in a first approximation (this technique is called dominant 
balance, Bender and Orszag 1978). Furthermore, for h(u) — ► 0, we expect w(n) to be more and more 
concentrated around fi = —1, so that in this range we can approximate the term D(fi,p)(l — /i 2 ) w 
2D_ 1 (l + /i), and I obtain w ps A 2 w with obvious solution tt; iwoexp (— A(/i + l)//i(u)). The factor 
X/h(u) can be determined by inserting this approximate expression for w into Eq. 3a. A trivial 
computation yields X/h(u) = 1/(1 — u). 

Now, going back to the equation for g, the spatial part of /, we find, also using the above, 

1 dg 2D_! 3 



, 2 - SyD-i (7) 
(/ az 7(1 — uj z 

from which, in the end, I find an approximate solution for the distribution function in the limit 
u 1: 

/ « ^lexp (8 7 3 J D_ 1 z) exp (-(/x + 1)/(1 - «)) . (8) 

It is thus seen that the detailed shape of the pitch angle scattering function D(n,p) is irrelevant, 
and that what is left of it (its value D_i at fi = — 1) only enters the spatial part of the distribution 
function /, not the angular one. 
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2.2. Downstream 

We make here the usual assumption, that the distribution function depends upon the particle 
momentum p as / oc p~ s in either frame (but see the Discussion for further comments). From 
the condition of continuity of the distribution function at the shock, denoting as p a and fj, a the 
particle's momentum and cosine of the pitch angle in the downstream frame, we have 

— w a (na) oc — exp(-(^t + 1)/(1 - u)) (9) 

Pa P 

where the irrelevant constant of proportionality does not depend on p,p a , fi, fi a . Using the Lorentz 
transformations to relate p,p a , H, Ma (/•* = (A*a ~~ u r)/(l — u r fi a ), p = p a 7 r (l — u r [x a ), with u r and j r 
the relative speed and corresponding Lorentz factor between the upstream and downstream fluids), 
I find 

w (a*o) = Tj Tjexp — -J. (10) 

For u — > 1, it is easy to derive from Taub's conditions (Landau and Lifshitz 1987) that u r — ► 1, and 
that (1 — u r )/(l — u) « 7 2 /7r - ► 2. This result does use a post-shock equation of state p = p/3, 
which is surely correct in the limit u — > 1. In the end, I obtain 

Wa{Ha) = , ~ , s exp (-2^-) . (11) 

(l-^ a ) s V l-Ato/ 
This equation shows why we needed to determine the pitch angle distribution, in the upstream 
frame, even for l-u/0: in fact, even though the angular distribution in the upstream frame (Eq. 
8) tends to a singularity, the downstream distribution does not (because the factor (1 — u)/(l — u r ) 
has a finite, non-zero limit), and the concrete form to which it tends depends upon the departures 
of the upstream distribution from a Dirac's delta. 

From now on I will drop the subscript a in /x a , since all quantities refer to downstream. In 
order to determine s, we now appeal to a necessary regularity condition which must be obeyed by 
the initial (i.e., for z = 0) pitch angle distribution, Eq. 11. Looking at Eq. 1 specialized to the 
downstream case, where u = 1/3 for very fast shocks, we see that this equation has a singularity 
at fj, = —1/3. Passing through this singularity will fix the index s. It is not convenient to use / 
directly; rather, I use its Laplace transform 

r+oo 

f(r,v,)= / f(z,»)e- rz dz. (12) 
J o 

Taking Laplace transforms of both sides of Eq. 1 I obtain 

- #t : KM + 7(l/3 + „)/= l± (d M(1 . (13) 

I am interested in the limit r — ► +oo. In fact, here I can use two results. First, in this limit, 
it is well-known (Watson's Lemma, Bender and Orszag 1978) that Eq. 12 reduces to 

/>, M )-» rt* = 'g> = ^. (14) 
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Despite this wonderful result in all its generality, I will actually use it only in the neighborhood of 
p = —1/3; here, Eq. 13 takes on a simple form: defining t = p, + 1/3, 

_ +at/ = -_ + -_ (15) 
r r ot 2 r at 

where I defined b = (jw a (p)D(p,p)(l - /J 2 ) ) | M= _i/ 3 , a = (jD(p,,p)(l - p 2 ))\^=-i/ 3 , and c = 

(d/dpD(p,p)(l — p?))/ D(p,p){l — /x 2 ) 1^= 1/3 - Now I make the Ansatz (to be checked a posteriori) 

that the term oc df/dt is negligible compared to the second order derivative in the limit r — > +oo. 
I am interested only in the most significant term in r, since Eq. 14 was only obtained to this order. 
Then, the equation 15 becomes 

atf = rW • (16) 
The above equation is the prototype of the one-turning point problem. Its solution, strictly in the 

neighborhood of the point t = p + 1/3 = 0, is (Bender and Orszag 1978, Sect. 10.4, Eq. 10.4.13b): 

f(t) « r^CAi^at) (17) 

where C is an arbitrary constant, and Ai(x) is one of Airy's functions. From this it can easily be 
checked that our Ansatz was justified. 

Clearly, close to the point i = // + 1/3 = 0, Eq. 14 and Eq. 17 must give the same results. 
Thus I find that, close to p, = —1/3, 

w a (p) oc Ai(r^ 3 a(ji + 1/3)) (18) 

which solves our problem: from this in fact we see that, since d 2 Ai(x) / dx 2 = xAi(x) by definition, 
and thus d 2 Ai(x)/dx 2 = in x = 0, then we must have 

d2yJa \ -n na\ 

^U=-i/3-0- (19) 

This is our sought for extra condition for s; we have seen that it comes directly from demanding 
that the boundary condition of the problem, Eq. 11, manages to pass through the singular point 
of Eq. 1, which I showed to be a conventional one-turning point problem familiar from elementary 
quantum mechanics. 



Q2„ 



dfi 2 



By substituting into Eq. 11 I find 

Wn -) =2- 2 ^3 2+s e- 1 (s 2 -5s + 3) = (20) 

/ ju=-l/3 

from which we obtain s = (5 ± \/l3)/2. The solution with the minus sign is unacceptable: in fact, 
if we plug Eq. 11 into the conservation equation 2, we see (Fig. 1) that for s < 3 the integral is 
< 0. We remarked after Eq. 3b, however, that this integral had necessarily to be > 0, so that we 
may conclude that 3 is an absolute lower limit to s. Thus we discard the solution with the minus 
sign, and are left with the unique solution 

s = — » 4.30 . (21) 




Fig. 1. — The integral of Eq. 3b with u = 1/3, for the angular distribution of Eq. 11, as a function 
of the parameter s, with arbitrary vertical scale. 



- 7- 



2.3. Oblique shocks 



Let us call 4> the angle that the magnetic field makes with the shock normal, in the upstream 
fluid. Then shocks can be classified as either subluminal or superluminal, depending upon whether, 
upstream, u/coscf) < 1 or u/coscf) > 1, respectively (de Hoffmann and Teller 1950). We are 
interested in the limit u — ► 1, so that most shocks will be of the superluminal type. In this case, we 
could (but we won't) move to a frame where the magnetic field is parallel to the shock surface, both 
upstream and downstream. However, downstream this extremely orderly field configuration appears 
more idealized than warranted by physical reality and observations. In fact, behind a relativistic 
shock, a large number of processes (compression, shearing, turbulent dynamo, Parker instability, 
two-stream instability) can generate magnetic fields; furthermore, there is no obvious reason why 
these fields should have large coherence lengths. In GRBs, a large number of observations of 
different afterglows supports this picture, the most detailed of all being those of GRB 970508, 
extending from a few hours to 400 d after the burst (Waxman, Frail and Kulkarni 1998; Frail, 
Waxman and Kulkarni 2000, and references therein). Accurate and successful modeling fixes the 
post-shock ratio of magnetic to non-magnetic energy densities to es ~ 0.1. Notice that here 
the protons' rest mass is not even the largest contribution to the non-magnetic energy density! 
Polarization measurements also support, albeit less cogently, the idea of a small coherence length: 
of the four bursts observed so far, only one has a detected polarization, at the 1.7% level (GRB 
990510, Covino et al, 1999). 

Thus the most plausible physical model downstream, is that particles move in a locally gener- 
ated turbulent, dynamically negligible magnetic field; if then we call I the average post-shock field 
coherence length, and restrict our attention to particles with sufficiently large energies (i.e., with 
gyroradii r g > I), we see there can be no reflection as particles approach the shock from down- 
stream. It follows that we expect the situation downstream to be identical to that of pure pitch 
angle scattering. Upstream, the parallel magnetic field is also irrelevant. In fact, backward deflec- 
tion of a particle occurs on a length-scale r g , but backward diffusion of the particle by magnetic 
irregularities only requires the sideways deflection by an angle ~ 7 -2 (7 being the shock Lorentz 
factor), for the shock to overrun the particle. This typically occurs on a length r]r g /j 2 , with rj a 
few. So, as 7 — ► 00, the length-scale for scattering upstream by the magnetic field increases, while 
that by magnetic irregularities decreases: the field is irrelevant. In the end, the same analysis as 
for pure pitch angle scattering applies, and the same index s and pitch angle distributions at the 
shock follow. 

In the case of subluminal shocks, a similar comment applies. Downstream, we expect on a 
physical basis the same situation as for superluminal shocks. Upstream, Eq. 1 is replaced by (Kirk 
and Heavens 1989) 



to which the same analysis as in Subsection 2.1 can be applied. Thus we find the same s and pitch 
angle distributions at the shock as above. 




(22) 
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As a corollary, it may be noticed that the above argument also implies that the results above 
are independent of the ratio the cross-field and parallel diffusion coefficients. 

3. Discussion 

For ultrarelativistic particles the energy spectral index k is related to s by 

l + \/l3 , N 

k = s - 2 = j — « 2.30 , (23) 

which is our final result. Also, now that we know s, the final pitch angle distribution at the 
shock, but downstream, can be determined (Eq. 11), and is plotted in Fig. 2. None of these 
results depends upon the specific form of D(fj,,p), so that widely differing assumptions should yield 
precisely the same results. 

How does this compare with numerical work? The near-constancy of the index s (or k) 
explains moderately well the results of previous authors: Kirk and Schneider (1987) find s = 4.3 
for their computation with the highest speed, which is however a modest Lorentz factor of 7 = 5, 
and a single functional dependence of D on /j,. Heavens and Drury (1988) find again a result of 
s = 4.2 for equally modest Lorentz factors, but for two different recipes for D. Extensive numerical 
computations using a MonteCarlo technique (i.e., totally independent of the validity of Eq. 1) 
were performed by Bednarz and Ostrowski (1998), for a wide variety of assumptions about the 
scattering properties of the fluids. They remarked quite explicitly that the energy spectral index 
k seemed to converge to a constant value, independent of shock Lorentz factor (provided 7> 1), 
magnetic field orientation angle (ft, and diffusion coefficient ratio, k±/k\\. They found k ~ 2.2. The 
present work confirms (for all untried forms for D(fi,p)) and extends their simulations (by yielding 
the exact value, and explicit forms for the particle angular distributions). The upstream angular 
distributions also agree well: Fig. 3a of Bednarz and Ostrowski clearly shows that this is (for the 
highest displayed value of 7 = 27) a Dirac's delta, in agreement with the largely limit in Section 
2.1. However, the downstream pitch angle distributions (in their Fig. 3b) agree well with mine, 
but not perfectly. Gallant, Achterberg and Kirk (1998) have claimed that there is a small error 
in Bednarz and Ostrowksi's distributions. As a matter of fact, my distribution (Fig. 2) agrees 
much better with Gallant et aVs and Kirk and Schneider's (1987) than Bednarz and Ostrowski's, 
despite the very small shock Lorentz factors of these two papers (7 = 2.3 and 7 = 5, respectively). 
Possibly, the small error in question may even explain the (small!) discrepancy between the two 
values of k. 

A limitation applies to the claim of universality of Eqs. 8, 11 and 23: I neglected any process 
altering the particles' energy during the scattering. Clearly, the results of this paper only apply in 
the limit dp/p < 1, where dp is the typical momentum transfer in each scattering event. In the large 
momentum limit considered here, it seems unlikely that this constraint may be violated. 

Lastly, a comment on the assumed dependence oc p~ s of the distribution function upon particle 




Fig. 2. — Pitch angle distribution at the shock in the downstream frame, Eq. 11 with s from Eq. 
21, with arbitrary vertical scale. 
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momenta is in order. It can be seen from Eq. 1 that such a dependence is not required by this 
equation. To see this, let us make the usual assumption that D is homogeneous of degree — r in p, 
i.e., D(fi,p) = q(p)p~ r - Then by defining a new variable z = z/p r , we see that the form assumed by 
Eq. 1 after this change of variable is identical to the original one, except that now p has altogether 
disappeared. At large z (i.e., far downstream), / — > = constant, and there is no p-dependence. 
This paradox is solved by noticing that the real problem to be solved involves both scattering (= 
Fermi acceleration) and injection. In this case, a typical injection momentum po arises naturally, 
and the dimensional problem discussed above is naturally solved: we must have / = f(...,p/po, ...) 
where the dots indicate all other parameters. In the limit of po — > 0, / does not tend to a constant 
independent of po as is always assumed, but tends instead to zero as / — > (po/p) s . Problems of this 
sort, though rare in astrophysics, are common in hydrodynamics, where they are called self-similar 
problems of the second kind (Zel'dovich 1956). They range from the deceptively simple laminar 
flow of an ideal fluid plast an infinite wedge (Landau and Lifshitz 1987) to the illuminating case 
of the filtration in an elasto-plastic porous medium (Barenblatt 1996). It is remarkable that, in 
the problem at hand, no such complication is necessary to fix the all-important index s, yet the 
powerful methods of intermediate asymptotics (Barenblatt 1996) and the renormalization group 
(Goldenfeld 1992) can be brought to bear on the intermediate 7 cases, where no easy limiting 
solution can be found. 

In short, what I have done in this paper is to show that the spectrum of non-thermal particles 
accelerated at relativistic shocks is universal, in the sense that the energy spectral index k, and 
the angular distributions in both the upstream and downstream frames (Eqs. 8, 11, 23, and 
Fig. 2) do not depend upon the scattering function D(fj,,p), the shock Lorentz factor (provided 
of course 7 > 1), the magnetic field geometry, and the ratio of cross-field to parallel diffusion 
coefficients. Thus we have the result that the cosmic rays' spectra are independent of flow details in 
both the Newtonian (Bell 1978) and the relativistic regimes. 
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